Fluctuation-driven insulator-to-metal transition in an external magnetic field 
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We consider a model for a metal-insulator transition of correlated electrons in an external magnetic 
field. We find a broad region in interaction and magnetic field where metallic and insulating (fully 
magnetized) solutions coexist and the system undergoes a first-order metal-insulator transition. A 
global instability of the magnetically saturated solution precedes the local ones and is caused by 
collective fluctuations due to poles in electron-hole vertex functions. 
PACS numbers: 71.30.+h,75.30.Kz,71.28.+d 



Metal-insulator together with magnetic transitions 
belong to the most important collective phenomena 
caused by the (screened) Coulomb electron repulsion. 
Particularly the correlation-driven Mott-Hubbard transi- 
tion from a paramagnetic metal to a paramagnetic insula- 
tor, i. e. with no apparent symmetry breaking, has been 
studied intensively in an effort to explain qualitatively 
the experimental findings in vanadium oxides doped with 
chromium Q . Although observed only at finite tempera- 
tures, the Mott-Hubbard transition is theoretically stud- 
ied mostly as an extrapolation to zero temperature and 
in the mean-field limit with only local dynamical fluctu- 
ations. Antifcrromagnetic long-range order is supposed 
to be suppressed at low temperatures by a frustration 
destroying the perfect-nesting property [§) . 

But even in this maximally simplified situation one is 
unable to solve the problem exactly. A scenario for the 
Mott-Hubbard transition in the spin-symmetric case was 
initially deduced from the so-called iterated perturbation 
theory (IPT) §, but it was later supported by analytical 
Q and numerical || renormalization-group arguments. 
The detailed transition behavior remains, however, still 
controversial ||. 

Unfortunately we do not have many cogent methods 
for the description of the critical behavior at the metal- 
insulator transition. The existing solutions are extrapo- 
lations. There is not a reliable interpolation between the 
Fermi-liquid regime and the strong-coupling, atomic like 
solution. IPT seems the only technique based on per- 
turbation theory being able to cover the desired features 
of the weak and strong coupling limits and to describe 
qualitatively the Mott-Hubbard transition. 

The description of the Mott-Hubbard transition in the 
spin-symmetric phase is further plagued with the nonex- 
istence of a parameter controlling the critical behavior. 
Nozieres recently suggested to extend the d — oo de- 
scription of the Mott-Hubbard transition to the magnetic 
case. The density of states (DOS) is no longer pinned at 
the Fermi energy and IPT does not reproduce correctly 
the strong-coupling limit ||. New methods are to be 
used to examine the magnetic case. 

On the other hand, magnetic field and spin polar- 
ization offer a natural parameter, magnetization, with 



the aid of which we can control the transition behavior. 
Moreover, the exact insulating fully spin-polarized state 
is explicitly known, namely the Hartree solution with well 
separated bands for up and down spins. The question 
whether the saturated magnet decays to a metallic state 
or whether the insulating phase remains stable for non- 
saturated magnetizations at intermediate and strong cou- 
pling can be studied as an extension of the Mott metal- 
insulator transition. 

A picture for the Mott transition near the fully satu- 
rated magnetic state in finite dimensions was suggested 
in || . It was shown that the fully saturated ferromagnet 
decays at B = B S (U) differently at weak and strong cou- 
pling. While in the former it is noncritical, a two-particle 
function gets critical at the ferromagnetic boundary in 
the latter. Recently this analysis was extended to the 
mean-field limit of a simplified Hubbard model with one 
spin species static |l(J. 

Unlike the spin-symmetric case, it is the metallic side 
of the transition that is difficult to describe in the mag- 
netic case. There is no complete solution below the ferro- 
magnetic boundary at intermediate and strong coupling. 
Only recently it was shown that for B /* B s the parquet 
approximation reproduces the exact asymptotics of the 
Bethe ansatz solution in d = 1 |IT |. However, a closed 
solution beyond weak-coupling Jl2[l3[ has not yet been 
constructed. 

The aim of this work is to investigate the exis- 
tence and kind of the fluctuation-driven insulator-to- 
metal transition in an applied magnetic field within the 
exact limit of high spatial dimensions. We concentrate 
on nearly ferromagnetic states at intermediate coupling. 

Using field-theoretic, diagrammatic techniques we 
show that at weak coupling, U < U Co s» 1.38w (for a 
semi-elliptic density of states with the bandwidth 2w), 
the insulating, fully magnetized solution goes continu- 
ously over to a metallic one by merging the bands of the 
up- and down-spin electrons at B S (U) = w — U/2. The 
behavior is controlled by one-particle functions. At inter- 
mediate and strong coupling, U > U CQ , however, the mag- 
netically saturated solution gets globally unstable along 
Bq(U) > w — U/2. A second critical point U Cl w 1.65 is 
found from which the saturated magnet gets also locally 
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unstable with a diverging electron-hole vertex function 
at B s > w - U/2 but B s < B a . 

The metallic solution at intermediate coupling is sta- 
bilized by correlations and the spectral function shows 
two satellite bands and a Kondo-like narrow quasiparticle 
peak around the Fermi energy. The metal remains stable 
also for B > B 8 (U) where already the fully spin polarized 
solution exists. It is destroyed at an upper critical filed 
B U (U) with a jump in magnetization to the fully mag- 
netically saturated state. We hence find a rather broad 
area of parameters U, B for which the metallic and insu- 
lating solutions coexist as local minima of the free-energy 
functional and within which a first- order metal-insulator 
transition takes place. We have a situation with mag- 
netic hysteresis and discontinuities in magnetization as 
known from metamagnetic transitions [pL 4|J 1 5| ] - 

Unlike weak coupling, the insulator-to-metal transi- 
tion at intermediate and strong coupling is controlled 
by two-particle vertex functions and their singularities 
(poles). No weak-coupling theory suppressing vertex 
corrections is able to describe this fluctuation-driven 
metal-insulator transition and the metallic solution with 
a narrow quasiparticle resonance correctly. Only self- 
consistent theories containing at least infinite series of 
multiple singlet and triplet electron-hole scatterings can 
offer an appropriate qualitative picture of the behavior of 
correlated electrons near the metal-insulator transition at 
intermediate coupling. 

We start with the Hubbard Hamiltonian in an exter- 
nal magnetic field B 



ff = ^(e(k)- M + a J B)4 



kcr^ko- 



(1) 



with a — ±1 denoting the spin direction and \x the chem- 
ical potential. We are interested in half filling at which 
[i = U/2 due to the electron-hole symmetry. The one- 
particle propagator can be represented as 



rr j B + 7J" m 



e(k)-£ CT (k,z) 



(2) 



where E CT is a dynamical correction to the Hartree static 
self-energy absorbed in the magnetization m. The dy- 
namical self-energy is determined from the Schwinger- 
Dyson equation of motion 



S CT (fc) 

u 



(3) 



' k'q 



T a _ a (k, k'; q)G a (k + q)G- a {k' + q)G- a {k') 



where r CT _ CT is the full two-particle vertex. We use a four- 
vector notation k = (k, iuj n ), q — (q, iv m ) for fermionic 
and bosonic variables, respectively. 

The vertex function T a - a has to be determined from 
a diagrammatic skeleton expansion. To allow for multiple 
electron-hole scatterings we must take into considerations 



at least the FLEX-type approximations with singlet and 
triplet electron- hole channels Jll],[l6|]. The electron-hole 
vertices are for the singlet and triplet channels, respec- 
tively: 



T a ^ a (k,k';q) = ^ + x 
/C a - a (k,k';q) 



U 



+ X a _ a (k-k>) 

u 



l-X n (q)X u (q) 
with dimensionless electron-hole bubbles 

Xaa'(q) = -^J2 G ^ G -'( k + ^ 



(4a) 
(4b) 

(4c) 



In this and the other FLEX approximations, the vertex 
at the two-particle scattering events remains unrenor- 
malized. These approximations are hence reliable only 
at weak or moderate interaction strengths. At strong 
coupling we have to pass to a higher level, the par- 
quet approximation with renormalized vertices at scat- 
tering events. At intermediate coupling, for U ~ 2u>, the 
FLEX and the two-channel parquet approximations de- 
liver qualitatively similar pictures in the spin-symmetric 
situation [ ft7| . 

The electron-hole bubble functions from 
X aa i (q) , are real and negative at the Fermi energy and 
the vertex functions get singular and negative if 



1+A ffff ,(q,0) < 



(5) 



Equality indicates the existence of a pole at the Fermi 
energy and negativity an instability of the solution, i. e. 
a local energy minimum turns into a local maximum. 

Bipartite lattices with perfect nesting and one elec- 
tron per lattice site get usually first unstable in the singlet 
channel at the antiferromagnetic point q = (it, 7r, ...). 
Frustration due to the inclusion of next-to-nearest or ran- 
dom hopping destroys this instability and a fully frus- 
trated model is expected to have no preference for a par- 
ticular low-temperature long-range order [^) . In the next 
steps we use this idealization of a fully frustrated sys- 
tem within the dynamical mean-field theory where the 
four-momenta are substituted by Matsubara frequencies. 
Further on, we use a semi-elliptic model density of states 
of the d = oo Bethe lattice po(E) = \Jw 2 — E 2 , with 
w = 1, approximating a three-dimensional DOS with a 
finite bandwidth and sharp edges. 

As a first step we investigate a local instability of 
the Hartree solution that can occur only in the singlet 
channel, since the triplet bubbles vanish in the saturated 
state. It is easy to find the explicit representation for the 
singlet electron-hole bubble 

X U (0) = -4U | A [1 - Bl] 3/2 + B m m| 
AU I —Vl-x 2 y/{2B m -x) 2 -l (6a) 

7T 
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with the equilibrium magnetization 

2 r 



arcsin (B m ) + B m y/l- B% n 



(6b) 



yJTi. The saturated magnetic solu- 



Here B m = B 

tion (m = 1) gets unstable if cither B + U/2 < 1 or 
X^(0) < — 1. The former condition follows from (6b), 
i. e. from the one-particle magnetization. The latter is 
the existence criterion for a pole in the vertex function 

The one-particle instability is realized at weak cou- 
pling for U < U Cl . The magnetization decreases and the 
DOS at the Fermi energy increases continuously upon 
decreasing the applied field B, Fig. [l]. 



the Fermi energy accompanied by a narrow, correlation- 
induced Kondo resonance. Upon increasing the magnetic 
field, the metallic solution remains locally stable and pre- 
serves the three-peak structure of the spectral function. 
The side bands get more pronounced and the central reso- 
nance narrower with a decrease in maximum, Fig. ^. The 
metallic solution gets locally unstable at an upper critical 
field B U (U) at which it again discontinuously goes over 
to the saturated magnetic insulator. The upper critical 
field is determined by a one-sided (B \ B u ) divergence 
in the longitudinal (local) magnetic susceptibility deter- 
mined by the triplet electron- hole vertex K,. The central 
Kondo resonance develops from a maximum to a sharp 
edge for the side bands with a jump in the Fermi-level 
occupation at the transition point. 
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FIG. 1. DOS of |-electrons at weak coupling for different 
applied fields. 

At a critical point U Cl ~ 1.65 both conditions (^) are 
fulfilled simultaneously and for U > U C1 the two-particle 
spin-flip instability ( |6a| ) at B S (U) > 1 — 17/2 causes the 
saturated magnet to get locally unstable and to decay to 
a magnetically non-saturated solution. 

We use the FLEX approximation with all three chan- 
nels for the quantitative description of the non-saturated 
solution at intermediate coupling. The self-energy from 
the most important electron-hole singlet channel has the 
explicit analytic representation 

= C/ /° — +-+) m ^f-^H, 



X a _ a (u)' - UJ + ) 
1 + X a - a (u' - UJ+) 



Im G- a (fi/ + )\ 



(7) 



where uj + = uj + i0 + . Analogous formulae hold for 
the other channels where only second-order contribution 
must correctly be subtracted. 

We find that the saturated solution at B = B S (U) 
for U > U Cl decays discontinuously, with a jump in the 
magnetization and energy, to a metallic solution with 
pronounced side bands and a sharp maximum in DOS at 
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FIG. 2. DOS of T-electrons at intermediate coupling for 
different applied fields. The local instability for the saturated 
solution lies at B « 0.135 

From the explicit analytic form of the self-energy con- 
tribution (Q) one immediately sees that infinite order 
summations are necessary to obtain an instability of the 
Hartree, fully spin-polarized solution for B + U/2 > 1 . 
In the fully polarized solution both Gj(w+) and 
are real for uj < 0, since the J.-DOS is unoccupied. How- 
ever, a pole in the vertex function T^i(uj+), i.e. when the 
condition 1 + A cr _ cr (ijj+) < is fulfilled for uj < 0, gener- 
ates a delta-function-like peak near the Fermi energy in 
the integrand. A nonvanishing contribution to E CT (i.e. 
the self-energy beyond Hartree) emerges even when start- 
ing the iterations from the polarized state. Clearly this 
would not occur in any weak-coupling treatment (e. g. 
second-order perturbation theory in U). Self-consistence 
is then needed to establish a locally stable solution. 

The coexistence region with two locally stable solu- 
tions, metallic and insulating, emerges due to a global 
instability of the magnetically saturated state. We find 
that instability condition ( pa|) is first fulfilled for a mag- 
netization m not obeying (pq) . The local two-particle 
instability at equilibrium is preceded by an instability 
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for a non-equilibrium magnetization. A first instabil- 
ity appears for m = 0, but for the triplet channel with 
Xf[. For the semi-elliptic DOS we then have an equation 
n/4U = §(1 - B 2 ) 3 / 2 defining an upper bound, B (U), 
for a global instability of the Hartree solution. The fully 
saturated magnetic state is hence globally stable only for 
B > B Q {U) and B > B S (U). 
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FIG. 3. Local (B 3 ) and global (-Bo) instabilities of the 
magnetically saturated solution with two critical points U co 
(■) and U C1 (•)• The dashed lines below the critical points 
show the respective instabilities of the non-saturated Hartree 
solution. 

However, the actual insulator-to-metal transition does 
not occur at the instability of either the insulating or the 
metallic solution. It is the energy that decides which 
solution is at global minimum. The line of the first- 
order metal- insulator transition, B c [u), depends on the 
approximation and lies in the coexistence region between 
the boundaries defined by local instabilities of the in- 
sulating and metallic solutions B s < B c < B u < Bq, 
respectively |18 . The line of first-order transitions ends 
up at a critical point U co where the line of the global, 
B {U), and of the local, B S (U), instability meet, Fig. ||. 

To conclude, we have investigated correlated electrons 
at zero temperature subjected to an external magnetic 
field and found that local instabilities of the insulat- 
ing, fully saturated magnetic solution at intermediate 
and strong coupling are preceded by a global instabil- 
ity. The saturated magnetic state decays via a first- 
order transition with a jump in magnetization to a non- 
saturated, metallic solution. We explicitly demonstrated 
the coexistence of metallic and insulating solutions on a 
mean-field (d = oo) model, but the same holds also for 
finite-dimensional systems near magnetic saturation. We 
have to take into consideration global fluctuations when 
investigating metal-insulator transitions and the stabil- 
ity and decay of fully polarized magnets. We further 
showed that the existence and shape of the metallic so- 
lution in the coexistence region is conditioned by the ex- 



istence of poles in two-particle electron-hole vertices. To 
understand the metal-to-insulator transition behavior in 
systems with correlated electrons completely requires to 
match the magnetic case studied here with the spin sym- 
metric situation and the Mott-Hubbard transition there. 
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